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Abstract. We have mapped the dust continuum emission from the molecular cloud covering a region of 28 pcx94 pc associated 

l/~) ■ with the well-known H n region RCW 106 at 1.2 mm using SIMBA on SEST. The observations, having an HPBW of 24" 

\Q ' (0.4 pc), reveal 95 clumps, of which about 50% have MSX associations and only 20% have IRAS associations. Owing to their 

higher sensitivity to colder dust and higher angular resolution the present observations identify new emission features and also 

show that most of the IRAS sources in this region consist of multiple dust emission peaks. The detected millimeter sources 

(MMS) include on one end the exotic MMS5 (associated with IRAS 16183-4958, one of the brightest infrared sources in our 

Galaxy) and the bright (and presumably cold) source MMS54, with no IRAS or MSX associations on the other end. Around 

O I 1 10% of the sources are associated with signposts of high mass star formation activity. Assuming a uniform dust temperature of 

20 K we estimate the total mass of the GMC associated with RCW 106 to be ~ 10 5 M Q . The constituent millimeter clumps cover 

a range of masses and radii between 40 to 10 4 M Q and 0.3 to 1.9 pc. Densities of the clumps range between (0.5-6) 10 4 cm -3 . 

We have decomposed the continuum emission into Gaussian and arbitrary shaped clumps using the two independent 
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CZ2 . 

structure analysis tools gaussclumps and clumpfind respectively. The clump mass spectrum was found to have an index a 

of 1.6 ± 0.3, independent of the decomposition algorithm used. The index of the mass spectrum for the mass and length scales 

covered here are consistent with results derived from large scale CO observations. 

X' 
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F~ i Key words. ISM: clouds - ISM: dust, extinction - ISM: H n and sub-millimeter wavelengths, large continuum maps essen- 

regions -ISM: structure - stars: formation tial for deriving a census of dust emission peaks arising from 

the deeply embedded phases of star formation are becoming 

available (Johnstone et al. 2000; Motte et al. 2001, 2003). 

1. Introduction RCW 106 was discovered by Rodgers, Campbell, & 

... , ,,. j. c , , .... , Whiteoak (I960) in a survey of Ha line emission of the south- 

Massive stars are believed to form in clusters within molecu- v ' J 

, , ■ , ,-. • . ., • • ,• . . ern Galactic plane. The Giant Molecular Cloud (GMC) asso- 

lar cloud complexes. Owing to their occurrence in more distant F v ' 

, , . „ , . , _. f .... , ■ -, ciated with RCW 106, was first detected by Gillespie et al. 

crowded stellar clusters, shorter formation timescales and for- ' J F 

... • f u- u ■ t t - t - »u e ..• e (1977) during observations of the molecular clouds associated 

mation in regions of high visual extinction, the formation of v ' & 

.„ , , . ■ r>. ill with southern Galactic H n regions in the J=l-0 transition of 

massive stars is still poorly understood. Dense, cool molecular & 

, , . ■ , • . . , , CO. Radio continuum mapping observations at 408 MHz and 

gas and dust cocoons enshroud massive protostars and account r± ^ & 

e . c ., .. .. . ■ ., • rp, ,. 5 GHz by Goss & Shaver (1970) detected a number of bright H 

for most of the extinction toward these regions. The extinction J v ' 6 

ii u- u <u t ■ t t »u a- »i n regions occurring along a line almost parallel to the Galactic 

is generally so high that a massive protostar cannot be directly & & & r 

, , . ., v u a /o o \ wi t t i /mm\ plane. IRAS observations have subsequently shown all of these 

observed even at the K band (2.2 //m). Watson et al. (1997) F M J 

, . ■ f . u j. . , .• f ., ■ ■ ■ H n regions to be bright infrared sources. The GMC is lo- 

succeeded in one of the rare direct observations of the ionizing & fe 

■ \r u a e t u ITPU • r, nn n£ a at r\ t - n cated at a distance of 3.6 kpc (Lockman 1979) covering an 

star in K band for the UCH n region G29.96-0.02. Optically F v ' ° 

... •„• . , , •„• . .• • • -, area of 70 xl5' between /=332.5 to 333.7 and appears as a 

thin millimeter and submilhmeter continuum emission from FF 

j . ■ ,. ., .. e e large complex of bright mid-infrared sources in the Midcourse 

dust cocoons shrouding the sites of massive star formation is fe F & 

, , , , c , , , , . ,-. • . ., Space Experiment (MSX) images between 8 and 21 urn (20" 

a valuable probe of molecular cloud structure. Owing to the F F v ' & r- \ 

. , , . • ., , . .. . , • . -,,. . resolution). The region contains one of the brightest infrared 

recent developments in the detection techniques at millimeter / & & 

H sources in our Galaxy (Becklin et al. 1973), IRAS 16183-4958, 

Send offprint requests to: B. Mookerjea, e-mail: bhaswati@phl.uni- which harbors the H n region G333.6-0.22. In addition to the 
koeln.de IRAS maps at the mid and far-infrared wavebands, far-infrared 
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(FIR) balloon-borne observations (1' angular resolution) of the 
dust continuum at 150 and 210 fim (Karnik et al. 2001) iden- 
tified 23 emission peaks with dust temperatures between 20 
and 40 K. Largescale CO(1-0) observations by Bronfman et al. 
(1989) with an angular resolution of 9' obtained a global view 
of the molecular gas distribution in this region. However most 
of the remaining spectroscopic observations toward this GMC 
are primarily pointed mode observations of the emissions of 
the lines of CO (Gillespie et al. 1977; Brand et al. 1984), CS 
(Gardner & Whiteoak 1978; Bronfman et al. 1996), H 2 CO 
(Gardner & Whiteoak 1984), NH 3 (Batchelor et al. 1977), H 2 0, 
OH and methanol masers (Braz & Epchtein 1983; Caswell 
1997; Walsh, Hyland, Robinson, & Burton 1997) toward pre- 
detected H n regions and IRAS sources. These observations 
suggest ongoing high-mass star formation activity and en- 
hanced density cores associated with most of the strong emis- 
sion peaks. This region will also be observed as a part of 
the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE), a Spitzer Legacy Science Program (Benjamin 
et al. 2003) at 3.6, 4.5, 5.8, and 8.0 //m using the Infrared Array 
Camera (IRAC). 

In this paper we present the first sensitive, large-scale con- 
tinuum survey of the RCW 106 region at millimeter wave- 
lengths. We have used the 1.2 mm observations to identify 
95 clumps, about half of which have MSX counterparts. We 
discuss association of these clumps with sources detected in 
MSX and IRAS surveys in particular. We have used methods 
based on clumpf ind (Williams & Blitz 1993; Williams et al. 
1994) and gaussclumps (Stutzki & Giisten 1990) to derive the 
clump mass distribution in this region. 



2. Observations and data reduction 

The 1 .2 mm continuum observations were carried out with the 
37 channel bolometer array SIMBA (SEST Imaging Bolometer 
Array) at the SEST (Swedish-ESO Submillimeter Telescope) 
on La Silla, Chile during the night from July 4 to 5, 2002. 
A mosaic covering approximately 1600"x5400" (a x 6) was 
observed using the fast scanning method. The mapped region 
extends over 28 pcx94 pc. Altogether 8 individual maps were 
used to construct the large mosaic and the total observing time 
was about 9.5 hours. The zenith opacity during the observation, 
measured using skydips, ranged between 0.209 and 0.385. For 
all but the two submaps at the top, two coverages were pos- 
sible, resulting in differences in the noise levels between the 
northernmost part of the map and the rest of it. The residual 
noise in the final co-added mosaic is 23 mJy/beam for most 
of the southern part of the mosaic and is around 50 mJy/beam 
for the northernmost part. Uranus was mapped for calibration 
purposes and a calibration accuracy of 15% was achieved. The 
HPBW for these observations was 24". 

All data were reduced and analyzed using MOPSI 1 accord- 
ing to the instructions of the SIMBA Observer's Handbook 
(2002) and also using methods described by Chini et al. (2003). 
Figure 1 shows the 1 .2 mm dust continuum emission from the 



1 MOPSI is a software package for infrared, millimeter and radio 
data reduction developed and constantly upgraded by R. Zylka 



region associated with RCW 106 as observed by SIMBA, to- 
gether with the associated IRAS point sources. 

3. Results 

To locate and measure the properties of dust emission peaks 
or "clumps" as we shall refer to them in this paper (see 
Section 3.4), we have used an automated procedure based 
on a two-dimensional variant of the clump-finding algorithm, 
clumpfind (Williams et al. 1994). The clumps so identified 
may have any arbitrary shape, so that the procedure gives re- 
sults similar to what an inspection by eye would have done but 
being automated, has far better accuracy. We further imposed 
the conditions that in order to be identified as a genuine sub- 
mm source, a clump must have (1) a deconvolved size larger 
than 50% of the HPBW (24") and (2) an S/N ratio greater than 
5. Table 1 lists the positions, peak intensities, integrated inten- 
sities, derived physical properties (size, mass and density) of 
the 95 clumps that have been detected at 1.2 mm. 

3.1. Overview of mm sources 

To obtain a first estimate of the nature of the sources detected 
at 1 .2 mm we have derived associations of these sources with 
existing MSX and IRAS point sources. We have associated a 
millimeter source with an MSX (or IRAS or FIR) source if 
the MSX (or IRAS or FIR) source lies within a radius of 40" 
(90" or 60 ") from the millimeter source. We have also iden- 
tified different signposts of massive star formation viz., H n 
regions, masers (H 2 0, OH, methanol etc) associated with the 
clumps. In addition, the MSX point sources for which the log- 
arithm of the ratio of the Band A (8 //m) to the Band E (21 //m) 
flux densities, log(S&/S2i) > 0.57, are identified as candidates 
for UCH ii regions (R. Simon, priv. communication). This is 
a direct extension of the Wood-Churchwell criterion (Wood & 
Churchwell et al. 1989) for identifying UCH n regions based 
on the the IRAS 12 and 25 //m flux densities. Figures 2, 3 and 
4 show overlays of the MSX (greyscale) 8 //m image with the 
SIMBA contours, with the mm-sources marked together with 
all MSX-UCH ii candidates in the region. We note that this 
method of deriving associations could be severely affected by 
the differences in the resolutions and sensitivities of the differ- 
ent observations under consideration. Following this method, 
we identify two main types of clumps: 

Type A : Clumps with an infrared (IRAS or MSX or FIR) 
association. Table 2 lists the 46 type A clumps together with 
the associated IRAS and MSX point sources, H n regions, 
MSX-UCHII candidates and H 2 0, methanol or OH masers. 
Owing to the higher angular resolution of the 1 .2 mm dataset 
presented here in contrast to the IRAS beamsizes, on several 
occasions a single IRAS source is associated with two or more 
clumps. All sources except five (MMS40, MMS52, MMS61, 
MMS63, MMS80) have MSX associations. Almost all clumps 
with MSX and IRAS associations are also identified as can- 
didates for MSX UCH n regions. Only 7 type A clumps are 
found to have H n regions associated with them and only 5 of 
them (MMS5, MMS9, MMS39/40, MMS68 and MMS84) have 
associated maser activities. MMS 54, having only a FIR associ- 
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Fig. l.The 1.2 mm mosaic 
of the star forming re- 
gions around RCW 106 ob- 
served with SIMBA. The 
typical ltr rms noise is 
23 mjy/beam. The contours 
correspond to 0.1, 0.2, 0.5, 
1, 2, 4, 10, 12, 15, 20, and 
40 Jy/beam. All associated 
IRAS point sources are also 
indicated in the plot.The 
mapped region extends over 
28 pcx94 pc. Also shown in 
the plot is the SIMBA beam. 



ation, has a methanol maser associated with it. Type A sources mass or size-wise no different from the Type A sources (see 
with associated H n regions are thus the more evolved, massive Section 3.3). 
star-forming clumps of this region. 



Type B : The remaining 49 clumps with no infrared asso- 
ciations. These are pure sub-mm sources, some of which are 
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Fig. 2. Contours of 1.2 mm dust continuum overlayed on the MSX 
8 fim image for the northernmost sub-region. The rms noise for the 
map is 50 mjy/beam and the contours are 0.2, 0.5, 1.0, 1.5, 2.5, 5.0 
Jy/beam and 10 to 40 Jy/beam in steps of 5. The millimeter clumps 
(+) and MSX-UCH n candidates (filled triangles) are marked. 
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Fig. 3. Same as in Figure 2 only for the middle part of the mapped 
region. The contours are at 0.1, 0.2, 0.5, 1, 1.5, 2.5, 5 and 8.5 Jy/beam. 



3.2. Mass 

Table 1 also presents the masses of the clumps (M c ] ump ) de- 
rived from their total integrated flux densities (F v ) using the 
following relationship: 

F V D 2 




KA (SO0O) 

Fig. 4. Same as in Figure 3 only for the southernmost part of the 
mapped region. The contours are at for 0.1, 0.2, 0.5, 1, 1.5, 2.5, 5 
10 and 12 Jy/beam. 



This equation is valid since dust emission at 1 .2 mm is op- 
tically thin. 

In equation 1, B v (Tdust) is the Planck function at the dust 
temperature Ti ast , k v is the dust mass opacity coefficient, D 
is the distance. In the absence of complementary continuum 
data, we generally assume a dust temperature of 20 K which 
is typically found in star forming cloud complexes (Johnstone 
et al. 2000; Motte et al. 2003, e.g.). However for the clumps 
with associated H n regions we have assumed T^st to be 40 K 
(Karnik et al. 2001). For the dust mass opacity coefficient we 
adopt the value of K23oghz = 0.005 cm 2 g~' Preibisch et al. 
(1993) assuming a gas-to-dust mass ratio of 100, and k v - 
K230GHz(y/230GHz)^ ', with/? = 2. Assuming a uniform dust tem- 
perature of 20 K, the total mass of the region is estimated to be 
10 5 M o . 

A few additional details about the calculation of the mass: 

(1) Given the large population of UC-H n candidates in the 
region, 1 .2 mm continuum emission observed with SIMBA is 
most likely to be contaminated by the free-free emission from 
the H n regions. However as mentioned earlier, no radio con- 
tinuum image of comparable angular resolution is available for 
the entire region. For the clumps with associated H n regions 
the contribution of the free-free emission from the gas never- 
theless has been estimated using the available 5 GHz flux den- 
sities (Haynes, Caswell, & Simons 1979) and a power law de- 
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, for optically thin free-free emission. 



The contribution varies between 4 and 36% (see Table 1). For 
the clumps for which the radio continuum flux densities were 
not available, despite their association with an H n region, we 
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Fig. 5. Plot of mass as a function of the radius (convolved with the 
beam) for all detected clumps. The dashed vertical line corresponds 
to the minimum size that a clump may have (resolution limit) and 
the other dashed line shows the minimum mass that a clump must 
have in order to be lcr above the mass detection limit. The sensitivity 
limit corresponds to MocR 2 . Also shown is the MocR 2 3 as found in CO 
observations (Heithausen et al. 1998), extrapolated to masses higher 
than 100 M . 



have assumed a constant 20% contribution to the flux densities 
at 1 .2 mm. 

(2) Two other sources of uncertainty in converting the ob- 
served flux density at 1.2 mm into gas mass are (a) the dust 
to gas mass ratio which typically ranges between 100 and 
150 and (b) the dust mass opacity coefficient, which ranges 
between 0.003 (diffuse ISM) up to 0.02 cm 2 g" 1 (Krugel & 
Siebenmorgen 1994) for dense protostellar environments, de- 
pending on the size and structure of the dust grains. 

3.3. Sizes & densities 

At a distance of 3.6 kpc the HPBW of 24" of SIMBA cor- 
responds to a diameter of 0.42 pc. The radii of the clumps 
identified by clumpfind range between the resolution limit 
to as large as 1.9 pc, with a substantial fraction having radii 
larger than 0.6 pc. Thus most of these clumps are of the size 
of star forming regions, rather than the innermost star forming 
cores typically of 0.01 to 0.1 pc size which have been dealt 
with in similar mm/sub-mm continuum studies of regions like 
p Ophiuchi, Orion B, W 43, IC 5146 etc. (Motte et al. 1998; 
Johnstone et al. 2000; Motte et al. 2003; Kramer et al. 2003). 

Figure 5 shows a log-log plot of the estimated mass as a 
function of the radii for all the 95 millimeter clumps in this re- 
gion. Since the clumps are not point sources, the completeness 
limit, measured in total flux or mass depends on both the range 
of object sizes and the surface brightness profiles. This holds 
for all surveys that attempt to measure the integrated flux from 
extended objects. Figure 5 shows that most of the more massive 
and larger clumps also belong to Type A and are either H 11 



regions or candidates for UC-H 11 regions. Figure 5 does not 
show any clear distinction between the two types of clumps. 
There appears to be a distinct mass-to-radius correlation for all 
the clumps, although part of it probably arises from the resolu- 
tion and sensitivity limits. The continuous line shows the mass 
to radius correlation (MocR 2 - 3 ) seen in CO clumps (Heithausen 
et al. 1998) for 0.01 < R/pc < 10 and 10~ 4 <M/M < 10 2 . For 
comparison, the Type A clumps show a MocR 28±03 correla- 
tion, whereas the Type B clumps exhibit a MocR 2 - 3±0 - 2 correla- 
tion. This implies that within the uncertainties the mass-radius 
relation for all types of clumps agree well with the result from 
Heithausen et al. (1998) as well as the sensitivity limit. 

The slope found for the dense, small prestellar objects stud- 
ied by e.g. Motte et al. (2001), is much shallower: M oc R iA 
though this relation may likewise be affected by the size- 
dependent detection thresholds. However, as pointed out by 
Motte et al. (2001) a linear relation is expected for self- 
gravitating isothermal condensations. 

Based on the calculated mass and derived sizes, in Table 1 
we also present the estimated densities of the clumps. We find 
that the clumps have moderate densities (~ 10 4 cm -3 ) spread 
over a rather narrow range. The densities are similar to those 
observed in the dust continuum study by Kerton et al. (2001) as 
well as the typical densities found in CO emission studies. The 
estimated densities also show that the sources under considera- 
tion here are quite different from the high density (~ 10 6 cm -3 ) 
condensations which are considered to be prestellar candidates 
(Motte et al. 1998). 



3.4. Nature of the mm sources 

Blitz (1993) proposed an operational categorization for struc- 
tures of different spatial scales seen in the ISM. In this scheme 
three types of structures are defined, viz., clouds, clumps and 
cores. Clouds are structures with masses > 10 4 M and cores 
are regions out of which single stars (or binary systems) form 
with masses < 10 M and sizes < 0.1 pc. Clumps correspond to 
structures with masses and sizes intermediate between clouds 
and cores, from which stellar clusters are formed. The millime- 
ter sources detected around RCW 106 have radii between 0.3 
and 1.9 pc and have masses between 40 and 10 4 M . Thus we 
refer to these sources as clumps and the entire region mapped 
with a mass of ~ 10 5 M as Giant Molecular Cloud (GMC). 

Given the limitations involved in associating the clumps 
with the radio observations detecting various signposts of star 
forming activities (viz., H n regions, masers), the evolution- 
ary stages of the clumps identified here are not accurate. The 
case of MMS54 is an apt example. While its fainter compan- 
ion MMS57 has an MSX counterpart, the stronger millimeter 
source MMS54 has no IRAS or MSX associations, but was de- 
tected both at 150 and 210 pm by Karnik et al. (2001) and has 
an associated methanol maser (Ellingsen et al. 1996, observed 
with an angular resolution of 4')- This definitely suggests that 
MMS54 shows early signs of massive star formation and de- 
mands further high angular resolution sub-mm dust continuum 
observations to characterize its nature. 
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Broadly, we propose the following evolutionary stages for 
the clumps: 

(1) Type A clumps with associated H n regions (and some 
of them with maser activities as well) are the most evolved high 
mass star forming regions. 

(2) Type A clumps which are candidates for MSX-UCH n 
regions (refer to Section 3.1 for the criterion) are high mass star 
forming regions, less evolved than those in (1). 

(3) Type A clumps which have only infrared associations 
and do not qualify as UCH n candidates are possible sites of 
low mass star formation at early stages of evolution. 

(4) The purely mm clumps (Type B) are cold condensa- 
tions, which given their densities are either prestellar or tran- 
sient features of the molecular cloud. 

Most of the existing spectroscopic observations are aimed 
at the sources (belonging to category (1) of Type A) that are 
bonafide massive star forming regions. Thus for the newly de- 
tected clumps there are no available spectroscopic observa- 
tions which may be used to estimate their gravitational stabil- 
ity and hence help to judge the potential of these clumps to be 
prestellar objects. 

4. Clump mass distribution 

Molecular cloud structure can be mapped via mm spectroscopy 
of molecular lines, continuum (sub)millimeter emission from 
dust or stellar NIR absorption by dust. The first method gives 
kinematical as well as spatial information and results in a three 
dimensional cube of data, whereas the latter two result in two 
dimensional datasets. We use the second method to derive in- 
formation about the structure of the molecular cloud around the 
RCW 106 region. 

For this purpose we have made use of the two most direct 
structure analysis tools, viz., clumpfind (Williams et al. 1994) 
and gaussclumps (Stutzki & Giisten 1990; Kramer et al. 
1 998) both of which decompose the observed emission into dis- 
crete clumps. The masses derived using clumpfind are already 
presented in Table 1. We have also used gaussclumps which 
follows a completely different clump identification and extrac- 
tion algorithm. In contrast to clumpfind, the gaussclumps 
algorithm iteratively deconvolves the observed emission into 
Gaussian shaped clumps. Though this algorithm was origi- 
nally written to decompose 3-dimensional molecular line data 
sets, it can also be applied to slightly adapted continuum data 
without modification of the code. 2 Motte et al. (2003) already 
used this method to analyze a continuum data set of W43. 

For both clumpfind and gaussclumps, we have selected 
only those clumps that have deconvolved sizes larger than 50% 
of the spatial resolution of the SIMBA data for further analysis. 
This discards the very small clumps which contribute only a 
tiny fraction to the total mass. 



2 We added two 2-D planes to the data set, mimicking two addi- 
tional but empty velocity planes bracketing the continuum data set. 
We checked that the formal values for the center velocities and the 
velocity widths of the clumps which the algorithm finds, stay con- 
stant, as expected. Thus, the algorithm effectively works only in the 
2-D center plane containing the continuum data. 



The masses of the identified clumps range between 40 and 
1O 4 M (Table 3). It is worth noting that particularly for the 
low mass clumps there is no one-to-one correspondence be- 
tween positions and masses of individual clumps found by 
gaussclumps and clumpfind. The gaussclumps algorithm 
decomposes a larger fraction of the total flux into clumps 
(Table 3) and thus finds a larger number of clumps. However 
we point out here that the purpose of this structure analysis and 
construction of the mass spectrum is to obtain an overview of 
the statistical properties of the clump ensemble. Thus, we have 
used all clumps regardless of their nature and/or evolutionary 
stages to build up the mass spectrum. 

Figure 6 shows the resulting distributions in a log-log plot 
of the number of clumps found per mass interval. We find 
that the clump mass distributions have similar shapes. The 
slope of the distribution flattens when going from the clumps 
of ~ 1O 4 M to clumps of ~ 1O 3 M . At smaller masses of 
~ 19OM corresponding to the 10cr mass at T^ util = 20 K, the 
distributions are almost flat. There is a sharp turnover of the 
spectra at still lower masses. No clumps are found beyond the 
minimum lcrmass. The turnover at masses of less than the 10cr 
mass has often been found in mass spectra created from CO 
surveys. Heithausen et al. (1998) showed that this turnover is 
not an intrinsic property of the clump distribution but rather due 
to the finite resolutions and the finite sensitivity. Observations 
at higher resolutions, show a smooth continuation of the mass 
spectra without any change in slope. A linear least squares fit 
to the clump distribution dN/dM oc M~ a for clumps with 
M > 10cr, results in linear correlation coefficients between 80 
and 90% and spectral indices a of between 1.7 and 1.6. Both 
spectra yield the same index within the standard deviation of 
0.3. 

As pointed out in Section 3.4, star forming cores are not 
detected by our survey. This is an important aspect to keep in 
mind while comparing the results of the structure analysis pre- 
sented here with the results of dust continuum and molecular 
line studies. 

Table 3 compares the properties of the clump ensemble ob- 
tained for RCW 106 with those of other dust continuum map- 
ping observations studying the fragmentation of the clouds. 
We restrict the table only to those studies which consider 
clumps that are more massive than stars and more extended 
than the prestellar cores considered e.g. by Motte et al. (1998). 
Interestingly, the spectral indices of the clump mass distribu- 
tions from these studies lie at ~ 1.6. 

This is in contrast to the recent dust continuum maps ob- 
tained at high resolution for the nearby clouds pOphiuchus, 
Serpens, and NGC 2068/71 (Motte et al. 1998; Testi & Sargent 
1998; Johnstone et al. 2000; Motte et al. 2001). The proper- 
ties of the resulting clump distributions have been compiled by 
Motte & Andre (2001). In these studies, clumps have masses 
down to 0.02 M Q , sizes less than 0.05 pc (< 10000 AU), and 
typical densities of 10 6 cm~ 3 . These starless and gravitation- 
ally bound clumps are therefore good candidates for being the 
progenitors of protostars. The above studies show a common 
spectral index of the clump mass distribution of ~ 2.3. There 
are only a few exceptions: the studies of Sandell et al. (2001) 
and Coppin et al. (2000) show indices of ~ 1 .5. 
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Fig. 6. The mass spectrum for the clumps identified using clumpfind 
(dashed histogram) and gaussclumps (drawn histogram) respec- 
tively. Error bars represent the standard deviation of a Poisson distri- 
bution "VAiv and are only drawn for one spectrum for clearness. The 
dashed vertical line indicates the lcr (=19M G ) and lOcr masses given 
by the average noise level of the map of 50 mjy/beam, the beam size 
(24", i.e. 0.42 pc), and assuming Tdust = 20 K. 



A spectral index of ~ 2.3 matches with the Salpeter IMF 
(Salpeter 1955) which is thought to hold for stars > 1M 
(Kroupa et al. 1993; Scalo 1998). There are indications inpOph 
that at still lower masses, the dust continuum mass spectrum 
flattens out, as does the stellar IMF. 

Interestingly, the mass spectrum found here for 
40 <M/M < 10 4 is similar within the uncertainties to 
the mass spectra derived from Galactic large-scale CO maps 
covering the mass range 10~ 4 <M/M Q < 10 4 (e.g. Heyer et al. 
2001; Simon et al. 2000; Kramer et al. 1998, and references 
therein) all of which agree with 1.5 < a < 1.8. In a recent 
study of CO in the Antennae galaxies (NGC 4038/4039), 
Wilson et al. (2003) found again similar index of 1.4 for a 
mass range from 5 10 6 to 9 10 8 M Q of the entire cloud complex 
using clumpfind. 

Thus the spectral index of the clump mass distribution 
found here is independent of the algorithm used and is consis- 
tent with both dust continuum and CO mapping studies which 
probe structures of densities ~ 10 3 - 10 4 ctrT 3 as we find in 
RCW 106. Spectral line maps at high spatial resolutions are 
needed to examine the stability of these clumps and find out 
whether they are gravitationally bound or are transient features 
of the turbulent clouds. 



5. Summary 

The 1.2 mm dust continuum emission from the GMC asso- 
ciated with the well-known H n region RCW 106 has been 
mapped at an angular resolution of 24" using SIMBA. In all, 95 
mm clumps have been identified, more than half of which are 



associated with pre-detected IR(MSX or IRAS) point sources. 
The region contains a number of dense massive star forming 
clumps which are recognized as candidates for UCH n regions 
and are bright in the far-infrared and millimeter wavelengths. 
There exist a large number of molecular condensations which 
are only bright in the millimeter wavelengths and have no IR 
associations. MMS54 has no MSX or IRAS associations but 
was detected in the FIR by Karnik et al. (2001) and has an 
associated methanol maser; it is a strong candidate for a star 
forming region at a very early stage of evolution. The contin- 
uum emission has been resolved in to Gaussian and arbitrary- 
shaped clumps using the two independent structure decompo- 
sition tools gaussclumps and clumpfind respectively. The 
masses of these clumps cover a wide range between 40 to 
10 4 M , suggesting that while the entire region is similar to a 
GMC, the sources detected are more like star forming regions 
(protoclusters) and not cores which would give rise to individ- 
ual stars. The distribution of the clump mass was found to have 
a spectral index a of 1.6 + 0.3, independent of the decomposi- 
tion algorithm used. The spectral index of their distributions 
match well with those observed in CO observations covering 
the same mass range. 

The higher angular resolution of this survey in contrast 
to the previous IRAS and far-infrared observations allows the 
identification of a number of clumps which are potentially 
interesting for further detailed continuum observations in the 
far-infrared and sub-mm as well as spectroscopic observations. 
Owing to the southern location of the region, ground-based ob- 
servations will have to wait for the upcoming telescopes like 
APEX and later ALMA. 
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Table 1. Observed and derived parameters of the millimeter sources detected in the RCW 106 complex. Dust masses have been calculated for 
Tdust = 20 K and an emissivity /3 = 2. Column (5) gives the clump radii (convolved with the beam). Estimated free-free contribution to the 
1.2 mm flux densities from the H n regions are ' 23%, 2 20%, 3 36%, 4 12%, 5 4%, 6 18%. The masses for these sources are calculated assuming 
a temperature of 40 K. 



Source 


RA(2000) 


Dec(2000) 


r peak 


Radius 


F v 


Mass 


nH 2 








mJy/beam 


pc 


Jy 


M 


10 4 cm" 3 


MMS1 


16:22:39.8 


-50:01:00 


828 


0.81 


2.5 


956 


1.74 


MMS2 


16:22:28.9 


-50:01:41 


585 


0.90 


2.9 


1097 


1.45 


MMS3 


16:22:41.6 


-50:04:36 


375 


0.60 


1.0 


376 


1.68 


MMS4 


16:21:49.3 


-50:05:43 


379 


0.63 


1.3 


500 


1.93 


MMS5 


16:22:10.1 


-50:06:06 


40013 1 


1.88 


129.5 


15936 


5.44 


MMS6 


16:22:25.9 


-50:06:04 


370 


0.89 


2.7 


1003 


1.37 


MMS7 


16:22:29.2 


-50:06:52 


430 


0.89 


2.2 


834 


1.14 


MMS8 


16:22:01.1 


-50:09:59 


502 


0.51 


0.7 


274 


2.00 


MMS9 


16:21:20.3 


-50:09:45 


4392 2 


0.83 


8.6 


1096 


5.44 


MMS10 


16:22:01.2 


-50:10:31 


475 


0.47 


0.8 


308 


2.87 


MMS11 


16:21:55.4 


-50:10:54 


356 


0.49 


0.5 


203 


1.67 


MMS12 


16:21:38.8 


-50:11:20 


258 


0.47 


0.4 


154 


1.43 


MMS13 


16:22:22.1 


-50:11:17 


775 


0.65 


2.5 


928 


3.26 


MMS14 


16:21:52.9 


-50:11:19 


273 


0.64 


0.8 


285 


1.05 


MMS15 


16:21:40.5 


-50:11:51 


546 


0.93 


2.8 


1052 


1.26 


MMS16 


16:22:04.6 


-50:11:58 


1015 


0.80 


3.2 


1187 


2.24 


MMS17 


16:21:08.0 


-50:15:14 


578 


0.32 


0.3 


109 


3.21 


MMS18 


16:21:32.3 


-50:18:49 


195 


0.51 


0.4 


139 


1.01 


MMS19 


16:21:44.9 


-50:18:48 


183 


0.41 


0.2 


83 


1.16 


MMS20 


16:21:46.6 


-50:19:03 


178 


0.24 


0.1 


30 


2.10 


MMS21 


16:21:49.9 


-50:19:03 


161 


0.55 


0.4 


158 


0.92 


MMS22 


16:21:50.8 


-50:19:27 


198 


0.72 


0.7 


259 


0.67 


MMS23 


16:21:28.4 


-50:23:12 


163 


0.47 


0.4 


139 


1.29 


MMS24 


16:21:40.9 


-50:23:11 


564 


0.62 


1.6 


590 


2.39 


MMS25 


16:21:39.3 


-50:23:53 


761 


0.84 


3.5 


1311 


2.14 


MMS26 


16:21:27.6 


-50:24:56 


1460 3 


1.25 


11.2 


1151 


1.34 


MMS27 


16:21:32.6 


-50:25:20 


3672 3 


0.96 


16.5 


1692 


4.34 


MMS28 


16:21:45.3 


-50:27:04 


375 


1.09 


3.1 


1183 


0.88 


MMS29 


16:21:32.7 


-50:27:12 


6502 


1.27 


24.7 


9281 


4.38 


MMS30 


16:21:42.9 


-50:28:15 


524 


0.96 


3.2 


1198 


1.31 


MMS31 


16:21:00.9 


-50:29:15 


308 


0.58 


0.8 


282 


1.40 


MMS32 


16:21:06.0 


-50:29:39 


319 


0.67 


0.9 


346 


1.11 


MMS33 


16:21:18.6 


-50:30:25 


1095 


0.90 


3.5 


1325 


1.76 


MMS34 


16:21:06.9 


-50:31:55 


282 


0.56 


0.7 


248 


1.36 


MMS35 


16:21:18.8 


-50:33:37 


190 


0.43 


0.3 


109 


1.32 


MMS36 


16:20:31.8 


-50:33:48 


225 


0.57 


0.6 


233 


1.22 


MMS37 


16:21:13.8 


-50:33:45 


451 


0.72 


1.7 


624 


1.62 


MMS38 


16:21:54.9 


-50:34:40 


283 


0.99 


1.6 


620 


0.62 


MMS39 


16:21:03.7 


-50:35:23 


8925 4 


1.19 


32.6 


4584 


6.18 


MMS40 


16:21:01.2 


-50:35:54 


4054 4 


1.14 


18.3 


2580 


3.95 


MMS41 


16:20:31.8 


-50:36:04 


151 


0.67 


0.7 


267 


0.86 


MMS42 


16:20:39.4 


-50:36:36 


804 


0.99 


4.5 


1675 


1.67 


MMS43 


16:20:38.5 


-50:37:41 


369 


0.46 


0.9 


331 


3.29 


MMS44 


16:20:40.2 


-50:37:56 


373 


0.71 


1.6 


593 


1.60 


MMS45 


16:21:12.1 


-50:38:02 


173 


0.40 


0.3 


98 


1.48 


MMS46 


16:20:30.1 


-50:38:12 


479 


0.87 


2.3 


849 


1.25 


MMS47 


16:20:49.5 


-50:38:43 


3322 


1.29 


12.1 


4530 


2.04 


MMS48 


16:20:28.4 


-50:39:24 


298 


0.74 


1.1 


406 


0.97 


MMS49 


16:21:15.6 


-50:39:29 


763 


0.74 


2.6 


965 


2.30 


MMS50 


16:21:13.9 


-50:39:45 


762 


1.11 


4.4 


1664 


1.18 


MMS51 


16:20:52.1 


-50:40:51 


348 


1.11 


2.7 


1029 


0.73 


MMS52 


16:20:29.4 


-50:41:00 


741 2 


0.93 


3.4 


431 


1.52 


MMS53 


16:20:36.9 


-50:41:00 


1268 2 


0.71 


3.5 


448 


3.55 


MMS54 


16:21:35.8 


-50:41:11 


2765 


1.14 


9.6 


3598 


2.35 


MMS55 


16:20:44.5 


-50:41:23 


173 


0.55 


0.5 


173 


1.00 
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Table 1. continued 
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Source 


RA(2000) 


Dec(2000) 


r peak 


Radius 


F v 


Mass 


nH 2 








mJy/beam 


pc 


Jy 


M 


10 4 cm" 3 


MMS56 


16:20:39.4 


-50:41:48 


968 2 


0.86 


4.3 


548 


2.45 


MMS57 


16:21:37.6 


-50:42:49 


359 


0.85 


1.8 


691 


1.09 


MMS58 


16:21:21.6 


-50:43:22 


326 


0.90 


1.5 


575 


0.76 


MMS59 


16:20:30.3 


-50:43:33 


326 


0.84 


1.8 


665 


1.08 


MMS60 


16:21:24.2 


-50:44:01 


153 


0.26 


0.1 


38 


2.09 


MMS61 


16:20:37.8 


-50:44:04 


893 


0.90 


4.2 


1563 


2.07 


MMS62 


16:20:54.7 


-50:44:11 


426 


0.93 


2.3 


860 


1.03 


MMS63 


16:20:28.6 


-50:44:44 


258 


0.87 


1.7 


642 


0.94 


MMS64 


16:20:48.8 


-50:44:59 


160 


0.60 


0.5 


203 


0.91 


MMS65 


16:20:32.0 


-50:49:56 


233 


0.69 


0.8 


293 


0.86 


MMS66 


16:20:35.4 


-50:51:17 


159 


0.67 


0.6 


237 


0.76 


MMS67 


16:19:30.4 


-50:52:32 


328 


0.63 


0.7 


278 


1.07 


MMS68 


16:20:11.9 


-50:53:17 


12460 5 


1.70 


36.1 


5548 


2.56 


MMS69 


16:20:32.1 


-50:53:24 


170 


0.72 


0.7 


255 


0.66 


MMS70 


16:20:16.1 


-50:56:38 


482 


0.78 


1.4 


541 


1.10 


MMS71 


16:20:06.8 


-50:57:02 


603 


0.94 


2.1 


789 


0.92 


MMS72 


16:20:02.6 


-50:58:30 


155 


0.92 


1.2 


436 


0.54 


MMS73 


16:20:06.8 


-51:00:06 


439 


0.62 


1.0 


391 


1.58 


MMS74 


16:19:47.4 


-51:00:15 


662 


0.73 


1.9 


702 


1.74 


MMS75 


16:20:47.6 


-51:00:27 


417 


0.63 


0.9 


342 


1.32 


MMS76 


16:20:03.5 


-51:00:38 


424 


0.52 


0.7 


259 


1.78 


MMS77 


16:19:41.4 


-51:01:11 


176 


0.77 


0.8 


297 


0.63 


MMS78 


16:19:52.5 


-51:01:35 


1905 


0.74 


5.0 


1867 


4.45 


MMS79 


16:19:49.1 


-51:02:24 


1331 


0.75 


3.9 


1476 


3.38 


MMS80 


16:19:55.8 


-51:02:23 


670 


0.61 


1.7 


639 


2.72 


MMS81 


16:19:24.4 


-51:02:33 


154 


0.66 


0.7 


248 


0.83 


MMS82 


16:19:25.3 


-51:02:56 


164 


0.44 


0.3 


120 


1.36 


MMS83 


16:20:15.4 


-51:03:25 


263 


0.86 


1.2 


470 


0.71 


MMS84 


16:19:38.9 


-51:03:28 


24996 


1.48 


21.9 


2871 


2.45 


MMS85 


16:20:02.7 


-51:03:26 


247 


0.33 


0.2 


79 


2.12 


MMS86 


16:19:16.7 


-51:03:45 


195 


0.55 


0.5 


195 


1.13 


MMS87 


16:19:22.8 


-51:03:53 


158 


0.76 


0.8 


319 


0.70 


MMS88 


16:19:09.2 


-51:03:53 


746 


0.86 


2.8 


1063 


1.61 


MMS89 


16:19:54.2 


-51:03:59 


529 


0.68 


1.3 


477 


1.47 


MMS90 


16:20:00.1 


-51:03:58 


300 


0.42 


0.5 


177 


2.31 


MMS91 


16:19:56.7 


-51:04:06 


258 


0.26 


0.2 


71 


3.90 


MMS92 


16:19:11.7 


-51:05:30 


159 


0.52 


0.4 


158 


1.09 


MMS93 


16:19:09.9 


-51:06:17 


315 


0.64 


0.8 


316 


1.16 


MMS94 


16:19:43.1 


-51:08:23 


155 


0.85 


0.9 


353 


0.56 


MMS95 


16:19:10.0 


-51:09:21 


224 


0.42 


0.3 


109 


1.42 
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Table 2. Type A : clumps with Infrared associations 



Source 


IRAS 


MSX 


Hn 


UCHn 
candidates 


Masers 


MMS2 




333.6896-00.2066 
333.6931-00.2023 




N 
N 




MMS4 




333.5684-00.1648 
333.5681-00.1608 




N 
N 




MMS5 


16183-4958 


333.6044-00.2165 
333.6046-00.2124 
333.6056-00.2044 


333.6-0.22 


Y 
Y 


H 2 0, OH 


MMS7 




333.6316-00.2562 




N 




MMS8 




333.5497-00.2455 




N 




MMS9 


16175-5002 


333.4678-00.1591 


333.46-0.16 


N 


CH 3 OH, OH 


MMS10 




333.5360-00.2565 




N 




MMS15 




333.4749-00.2362 




Y 




MMS16 


16182-5005 


333.5328-00.2670 




N 




MMS17 




333.3755-00.2007 




Y 




MMS26 


16177-5018 


333.3080-00.3655 


333.3-0.4 


Y 




MMS27 


16177-5018 


333.3080-00.3655 


333.3-0.4 


Y 




MMS28 




333.3113-00.4254 




N 




MMS29 




333.2898-00.3898 




N 




MMS33 


16174-5022 


333.2271-00.4053 
333.2276-00.4084 




Y 

N 




MMS36 




333.0929-00.3539 




N 




MMS39 


16172-5028 


333.1313-00.4268 


333.1-0.4 


Y 


CH3OH 


MMS40 


16172-5028 




333.1-0.4 




CH 3 OH 


MMS45 




333.1241-00.4756 




N 




MMS47 




333.0688-00.4445 




N 




MMS49 




333.1082-00.5011 




Y 




MMS50 




333.1082-00.5011 




Y 




MMS52 


16167-5034 




333.0-0.4 






MMS53 


16167-5034 


333.0149-00.4431 


333.0-0.4 


Y 




MMS54 


S19* 






N 


CH 3 OH 


MMS56 


16167-5034 


333.0125-00.4555 


333.0-0.4 


Y 




MMS57 




333.1166-00.5816 




N 




MMS59 


16167-5036 










MMS60 




333.0675-00.5855 




N 




MMS61 


16167-5036 


332.9881-00.4864 




Y 




MMS63 


16167-5036 










MMS67 


16156-5046 










MMS68 


16164-5046 


332.8269-00.5489 


332.8-0.6 


Y 


H 2 


MMS70 




332.8007-00.5953 




Y 




MMS74 




332.7024-00.5866 




N 




MMS75 


16170-5053 


332.8137-00.6980 




Y 




MMS78 




332.6926-00.6121 




N 




MMS81 




332.6213-00.5701 




N 




MMS82 




332.6213-00.5701 




N 




MMS84 


16158-5055 


332.6555-00.6129 


332.7-0.6 


Y 


CH 3 OH,H 2 


MMS85 




332.6853-00.6463 




N 




MMS90 




332.6853-00.6463 




N 




MMS86 


16153-5056 










MMS88 


16153-5056 


332.5860-00.5634 




Y 




MMS92 


16152-5059 










MMS93 


16152-5059 


332.5600-00.5902 




N 





FIR detection at 150 and 210 yum by Karnik et al. (2001). 
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Table 3. Comparison of the masses of the clumps identified in RCW 106 with large scale dust continuum studies of molecular clouds finding 
clumps or fragments more massive than 30 M and sizes of more than 0.05 pc. N cl s the number of clumps identified. Sizes refer to deconvolved 
FWHM sizes. M t ( [ lmps is the total mass of all clumps, a is the spectral index of the mass distribution shown in Figure 6 and n is the density in 



Cloud 


Distance 


HPBW 


N ci 


Sizes 


Mass Range 


..clumps 
JM tot 


a 


log(n) 


Reference 




[pc] 


["/pc] 




[pc] 


[Mo] 


M Q 




log[cirT 3 ] 




RCW 106 


3600 


24/0.4 


95 


0.2-3.7 


37-16000 


1.0 10 5 


1.7 ±0.3 


3.7-4.5 


this paper using clumpfind 


RCW 106 


3600 


24/0.4 


122 


0.2-3.7 


110-14000 


1.5 10 5 


1.5 ±0.3 




this paper using gaussclumps 


KR140 


2300 


14.5/0.16 


22 


0.2-0.7 


0.5-130 


590 


1.5 ±0.04 


3.1-5.1 


Kerton et al. (2001) 


M8 


1700 


12/~ 0.1 


75 




0.3-43 


280 


1.7 ±0.6 




Tothill et al. (2002) 



